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Abstract. The use of Type la supernovae as calibrated standard candles is one of the most powerful tools to 
study the expansion history of the universe and thereby its energy components. While the analysis of some ~50 
supernovae at redshifts around z ~0.5 have provided strong evidence for an energy component with negative 
pressure, "dark energy", more data is needed to enable an accurate estimate of the amount and nature of this 



energy. This might be accomplished by a dedicated space telescope, the SuperNova / Acceleration Probe (200C; 
SNAP), which aims at collecting a large number of supernovae with z < 2. In this paper we assess the ability 
of the SNAP mission to determine various properties of the "dark energy." To exemplify, we expect SNAP, if 
operated for three years to study Type la supernovae, to be able to determine the parameters in a linear equation 
of state w(z) = u>o + Wi z to within a statistical uncertainty of ±0.04 for u>o and 1° *7 f° r Wl assuming that the 
universe is known to be flat and an independent high precision (an m ~ 0.015) measurement of the mass density 
Q m , is used to constrain the fit. An additional improvement can be obtained if a large number of low-z, as well 
as high-z, supernovae are included in the sample. 
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1. Introduction 

The description of the universe lies at the heart of cos- 
mology, and it is not surprising that several methods aim- 
ing at the determination of cosmological parameters cur- 
rently are considered. For example, the power spectrum 
of the cosmic microwave background radiation provides 
means to determine the total energy content of the uni- 
verse, for which recent results of the balloon-based CMB 



measurements (Jaffe et al. 200C) quote the value fitot = 
1.05 ± 0.040. Constraints on the matter energy density of 
the universe, f2 m , can be derived, e.g., from galaxy clus- 
ter abundances (Bahcall & Fan 1998| , Carlberg et al 1998Q , 
and large-scale structure (Peacock et al 2001). These tests 
are consistent with Q m ~ 0.3, see however (Blanchard et 
al. [2000 ). Furthermore, studies of weak lensing effects of 
background objects in mappings of the sky provides infor- 
mation about the mass distribution in the universe, and 
thus niGcisiirGS ^Ifji . Sgg, e.g., van Waerbeke et al. ( |1999[ ) 
for a discussion of the accuracy of this method. 
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1 This value was derived assuming that the Hubble constant 
is 71 ± 8 km s _1 Mpc _1 



On top of this, measurements of supernovae at var- 
ious redshifts provide a simple way to estimate cosmo- 



logical parameters (Goobar & Perlmutter 1995). In fact, 
this is the aim of at least two collaborations (Riess et 
al. 1998, Perlmutter et al. 1999), both of which recently 



have published data in favour of a large energy compo- 
nent attributable to a cosmological constant, or an evolv- 
ing scalar field such as "quintessence" (Ratra & Peebles 
1988| , Caldwell et al. |1998| ). The feasibility to determine 



the properties of this "dark energy" component by using 
supernova data has recently been considered by several au- 
thors (see, e.g., Huterer & Turner |l999] , Saini et al. 2000 



Maor et al. 2001, Astier 2001, Weller & Albrecht 2000 and 



Barger & Marfatia 200C , just to list a few), and concl usion s 
vary significantl y. For instance, Huterer & Turner ( 1999 ), 
and Saini et al. (200C) devise methods for reconstructing 
the potential of an acceleration-driving scalar field, using 
supernova measurements. On the other hand, Maor et al. 
( 2001 ) assess the possibility to use supernovae to distin- 
guish between various cosmological models, allowing for 
an evolving equation of state w{z) (which is equivalent to 
scalar- field models) . They conclude that the prospects for 
determining the equation of state in this way are bleak. 
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Barger & Marfatia (200C) support this latter view, exem- 
plifying how particular data realisations may give mislead- 
ing conclusions regarding the dark energy. Again, Weller 
& Albrecht (200C) are more optimistic regarding a de- 
termination of w(z), provided that accurate independent 
estimates of the matter energy density f2 m are at hand. 



I = 



As already emphasized by one of us (Astier 2001), much 
of the discrepancies stem from differences in the initial 
assumptions, e.g., in the prior knowledge of il m . 

In this paper we intend to study the extent to which 
properties of the dark energy can be determined, assum- 
ing that observations of a large number of supernovae at 
high redshifts become available. Such data could be pro- 
vided by the projected SNAP satellite mission. In section 
^ we establish our notation and give the expression for 
the luminosity distance dz. Section || contains investiga- 
tions of different scenarios in line with the SNAP pro- 
posal (200C). Confidence regions for cosmological parame- 
ters are obtained for various situations. Section |] considers 
the relative importance of events at various redshifts by 
investigating the effect of adding a small sample at vari- 
ous specific redshifts. In section [|, we analyse the system- 
atic errors in cosmological parameter estimation that are 
caused by gravitational lensing. We end with a discussion 
of the main conclusions in section ^. Appendix [A] outlines 
the construction of our log-likelihood functions in some 
detail. 

2. Apparent magnitude and luminosity distance 

We intend to investigate the feasibility to determine cos- 
mological parameters = (fl mi flx,wo, w i) by using ob- 
servational data from supernovae at different redshifts z. 
Here, f2 m and fix denote the present-day energy den- 
sity parameters of ordinary matter fi m (z) and a "dark 
energy" component fix(z), respectively. The equation of 
state w(z) of the dark energy is parametrised by {wq, w\) 
to linear order: w(z) « wq + w\ z. 

The apparent magnitude m of a supernova at redshift 
z, assuming the cosmology 9, is given by 



m{0,M,z) = M + 51og 10 [d' L {d,z)}, 
M = 25 + M + 51og 10 (c/iJ ), 



(1) 
(2) 



where M is the absolute magnitude of the supernova, and 
d' L = H{)dL is the HQ-independent luminosity distance, 
where Ho is the Hubble parameter 2 . Hence, the inter- 
cept A4 contains the "unwanted" parameters M and Hq 
that apply equally to all magnitude measurements (we do 
not consider evolutionary effects M = M(z)). The Hq- 
independent luminosity distance d' L is given by 



d'r 



(i + ^7=n^sin(7=n^/), n fe <o 

(1 + 2;)/, r2 fc = o 

(i + z)^ S inh(vn^/), n k >o 



(3) 



Qv. = 1 — \lrn, — ^ 



(4) 



In the expression for M, the units of c and Ho are km s 
and km s _1 Mpc^" 1 , respectively. 



dz' 



Jo H'(z'Y 
H'(z) = H(z)/H = 



f(z) = cxp 



y/(l + zY fi m + f(z) Q x + {1 + zf Q k , 
, l + w(z')~ 



3 / dz' 

o 



1 



where we consider an equation of state linear in z: 
w(z) = Wq + Wi z. 



(5) 

(6) 
(7) 

(8) 



3. Statistical uncertainties for one year of SNAP 
data 



The SuperNova/ Acceleration Probe fl2000[ SNAP) is a 
proposed two-meter satellite telescope specifically de- 
signed to discover and follow supernovae over a wide red- 
shift range. In particular, such an instrument would be 
able to provide photometry and spectra of more than 2000 
SN la per year (SNAP proposal 2000). We will investi- 
gate the accuracy of cosmological parameter estimations 
based on one year of SNAP data. To this end, we as- 
sume that 2000 supernovae are obtained in the redshift 
interval z £ [0,1.2], and an additional 100 at high red- 
shift z £ [1.2, 1.7]. The individual measurement precision 
is assumed to be Am = 0.15 magnitudes, including the 
intrinsic spread of supernova brightnesses. We divide the 
redshift interval into bins of equal size Az = 0.05. In sum- 
mary: 

z range # SNe # bins # SNe/bin prec./bin 

[mag.] 

[0.0,1.2] 2000 24 83.33 0.0164 

[1.2,1.7] 100 10 10 0.0474 

We will use the fiducial cosmology from the SNAP pro- 
posal: 0t. r uc = (0.28,0.72,-1,0). These assumptions ad- 
here to the SNAP proposal ( |2000 ), except that we do not 
include any systematic errors. However, in section || we 
will investigate the effects of gravitational lensing on cos- 
mological parameter estimations. 

Below, we consider several different scenarios and 
present confidence regions for parameter estimates. The 
methodology that has been employed is outlined in App. 

The one-parameter one-sigma un certa i nties for the var- 
ious cases are summarised in tables A.l - A. 3. 



3.1. Confidence regions for (O m , tlx) 

First, let us assume that it is known that the dark energy 
corresponds to a cosmological constant, so that (wq, w±) = 
(—1,0). In this particular case, Six is often denoted f^A- 
Figure |l] shows confidence regions for (VL m ,Slx) for vari- 
ous situations. As regards O m , we assume either no prior 
knowledge, or else prior knowledge with Q m Gaussian 
around the true value with asi^-prior = 0.05. Concerning 
the intercept M. , we assume either exact knowledge of M. , 
or no prior knowledge at all. The latt er ca se involves the 
expression y? M _- lvA , given in appendix A.l . 
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Fig. 1. 68.3 % confidence regions for (fl m ,flx) in the 
one-year SNAP scenario. The filled region assumes exact 
knowledge of M. (solid and dashed lines approximately co- 
incide). A full three-parameter fit with no prior knowledge 
of M. is assumed for the two larger confidence regions: 
the region with a dotted line assumes no prior knowl- 
edge of fl m , while the dash-dotted line assumes a prior 
knowledge with fl m Gaussian around the true value and 
O"n m -prior = 0.05. 



Fig. 2. 68.3 % and 95% confidence regions for (fl m ,flx) 
in the one-year SNAP scenario without the 100 events 
for which z £ [1.2,1.7]. The filled region (solid line) as- 
sumes exact knowledge of M. , and the dashed line within 
the filled region assumes also a prior knowledge with fl m 
Gaussian around the true value and en —prior = 0.05. 
A full three-parameter fit with no prior knowledge of M. 
is assumed for the two larger confidence regions: the re- 
gion with a dotted line assumes no prior knowledge of fl m , 
while the dash-dotted line assumes a prior knowledge with 
fl m Gaussian around the true value and cm m -r>rioi — 0.05. 



Under the assumption of exact knowledge of A4, we 
find the uncertainties in fl m and fix to be <Jn m ~ 0.015, 
(Tsj x 0.027. However, with no prior knowledge of Ai, 
the uncertainty in fix grows almost by a factor of two. 
Note that the uncertainty in fl m is essentially unaffected. 
Hence, imposing the prior knowledge of fl m as outlined 
above, does not significantly affect the size of the confi- 
dence region. To emphasize the importance of obtaining 
at least a few supernovae at high redshift, we perform 
the same calculation including only the events for which 
z < 1.2, see figure ||. Even though there were only 100 
such supernovae in the original calculation, they result in 
about 25 % better determination of fix- Thus, it seems to 
be well-worth the effort to devise a scheme for obtaining 
these high- z events. On the other hand, in order to reduce 
the sensitivity to uncertainty in the intercept , it is im- 
portant to have supernovae at low redshifts. To illustrate 
this, we have examined a situation where the redshifts of 
the 2000 supernovae at z £ [0, 1.2] are distributed accord- 
ing to a constant rate per co-moving volume element (as 
opposed to the uniform distribution used before). The few 
events for z £ [1.2,1.7] are still considered to be uniformly 
distributed. As seen in figure ||, this does hardly affect the 
uncertainties when M. is exactly known. However, for the 
worst-case scenario of no prior knowledge of M, the un- 
certainty in fix grows almost by a factor of three. The 
relative importance of events at various redshifts is fur- 
ther discussed in section H. 



3.2. Confidence regions for (fl m ,wo) or (Qx,wo) 

Next, we assume that the equation of state of the dark 
energy can be described by a constant w — wq, so that 
wi = 0. Figures f|- [|show confidence regions for (fl m , wq) 
or (fix, wo) under different assumptions that fix one pa- 
rameter in the expression f2 t ot = flm + fix ■ figure |] fixes 
the total energy density f2 to t = 1, which corresponds to a 
flat universe; figure |^ assumes that the density of the dark 
energy is known exactly, fix = 0.72; figure || assumes that 
the energy density of ordinary matter is exactly known, 
fl m = 0.28. As before, we consider either exact knowl- 
edge of M., or no prior knowledge. We also consider prior 
knowledge of either fl m or O to t) with spread <7 P rior = 0.05. 
It turns out that the (unrealistic) case where fix is well- 
known gives the best determination of the other param- 
eters under consideration, and that a good knowledge of 
^tot is preferred over a well-determined fl m . 

In figure fj], all three parameters (fl m ,flx,WQ) are al- 
lowed to vary, while both fl m and r2 to t are independently 
subject to Gaussian priors with <Tr2 m _p r i r = 

0.05 and 

°~Qtat— prior = 0.05. Comparing the case with exact M to 
the situation with no prior knowledge, it can be noted that 
the uncertainty of the latter mainly grows in wq. 



3.3. Confidence regions for (wq,wi) 

Recently, Maor et al. fl200l| MBS) considered the prob- 
lem of determining the equation of state of the dark en- 
ergy using supernova measurements. In particular, they 
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Fig. 3. 68.3 % and 95 % confidence regions for (f2 m , fix) 
in the one-year SNAP scenario with a constant rate per 
co-moving volume for z £ [0, 1.2], and the 100 z £ [1.2, 1.7] 
events uniformly distributed. The filled region assumes 
exact knowledge of M. (solid and dashed lines approxi- 
mately coincide). A full three-parameter fit with no prior 
knowledge of M. is assumed for the two larger confidence 
regions: the region with a dotted line assumes no prior 
knowledge of fi m , while the dash-dotted line assumes a 
prior knowledge with Q m Gaussian around the true value 
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Fig. 4. 68.3 % confidence regions for (O m , wo) in the one- 
year SNAP scenario with the flatness assumption f2 tot = 
1. The filled region (solid line) assumes exact knowledge of 
M. , the dashed line within the filled region assumes also a 
prior knowledge with fl m Gaussian around the true value 



and an„ 



-prior 



= 0.05. The dotted and dash-dotted lines 



assume no prior knowledge of A4, without and with fi r , 
prior, respectively. 



investigated an idealised experiment with thousands of 
supernovae in the redshift range z £ [0,2], divided into 
50 bins. The relative precision of the luminosity-distance 
d L was taken to be 0.6 % per bin, which corresponds to 
a magnitude precision <7i = 0.006 x 5/ln(10) ~ 0.0130 
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Fig. 5. 68.3 % confidence regions for (fl m ,Wo) in the one- 
year SNAP scenario assuming exact knowledge of fix, 
i.e., no prior knowledge on the geometry. The filled re- 
gion (solid line) assumes exact knowledge of M., the 
dashed line within the filled region assumes also a prior 
knowledge with Gaussian around the true value and 
&n m -prior = 0.05. The dotted and dash-dotted lines as- 
sume no prior knowledge of M, without and with fl m 
prior, respectively. 
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Fig. 6. 68.3 % confidence regions for (fix, Wq) in the one- 
year SNAP scenario assuming exact knowledge of fl m ■ The 
filled region (solid line) assumes exact knowledge of A4, 
the dashed line within the filled region assumes also a 
prior knowledge with f2 tot Gaussian around the true value 
and o"Q tot _p r ior = 0.05. The dotted and dash-dotted lines 
assume no prior knowledge of M, without and with f2 tot 
prior, respectively. 



for each bin. The equation of state is taken to be linear, 
w(z) = wq + wi z. Confidence regions for (u>o, w\) were de- 
termined, using the cosmology 8 tvnc = (0.3,0.7,-0.7,0). 
The log-likelihood was determined both for an exact f2 m , 
and with fi TO integrated over f2 TOj truc ± 0.1. 

Figure || shows our calculation of the confidence re- 
gions for this scenario. This figure should be compared 
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Fig. 7. 68.3 % and 95 % confidence regions for (Qx,wo) 
in the one-year SNAP scenario assuming prior knowledge 
with U, m and Qtot independently Gaussian around their 



true values and co 



-prior 



0.05 



en 



tot -prior ■ 



The hrst 



figure assumes exact knowledge of Ai, while the second 
assumes no prior knowledge of M. . 



with figure 2 of MBS. (Since MBS present one- and two- 
sigma contours, rather than the 68.3 % and 95 % confi- 
dence regions, we have included both cases to facilitate 
comparison.) There is a considerable discrepancy between 
these figures and MBS0. 

In conclusion, it seems to us that this scenario en- 
ables a better constraining of (wo,Wi) than was pre- 
viously anticipated by MBS. However, the scenario as- 
sumes that more than 6000 supernovae uniformly dis- 
tributed over a rather optimistic redshift range are ob- 
served. Consequently, in this section, we calculate (wq, w%) 
confidence regions for the cosmology of MBS, using the 
weaker precision and a smaller redshift range assumed in 
the SNAP proposal (2000), see figure p"2| below. However, 
we focus our attention on the fiducial cosmology of SNAP: 
truc = (0.28, 0.72, -1, 0) (see figures | - [ll] and |l|) . 

We consider the ability to determine the equation of 
state of the dark energy to linear order, w{z) =»o+wi z. 
We will assume flatness, J7tot = 1, and impose some prior 
knowledge of O m . Figure || shows confidence regions for 
various assumptions regarding M. and f2 m . We mainly 
consider a Gaussian prior with crf^ m _p r j or = 0.05. The uni- 
form prior with Q m £ f2 mjtr ue±0.1 is considered in the case 
of exact knowledge of M.. since this is the situation con- 
sidered by MBS. In figure [h], the few high-z supernovae 
have been excluded. When M. is exactly known, these are 
not so important in determining (wo ,w±) as they are for 
(fi m ,fix)i basically because fix becomes less significant 
with increasing redshift. However, note that the high-z 
events make some difference when M. is poorly known. 
Figure [ll] shows the situation when the supernovae at 
z G [0,1.2] are distributed according to a constant rate 
per co-moving volume element. As can be expected, un- 
certainties are not affected when M. is considered to be 
exactly known, but degrade considerably with no prior 
information of Ai . Figure |l2| considers the same scenario 



3 It has come to our attention that MBS used <Ji = 0.03 mag., 
which corresponds to a relative precision in cLl of about 1.4 %, 
and that their contours really correspond to 68.3 % and 95 
% confidence regions (Brustein, private communication). This 
fully accounts for the discrepacy between figures. 



as in figure ^ as regards precision and priors for fl m , but 
uses the fiducial cosmology of MBS. 

With the priors for Q m assumed in figures |9| - [l2], 
the equation-of-state parameters are rather poorly con- 
strained by one year of SNAP data, especially when M. is 
left unspecified. In order to see what SNAP can achieve 
over its expected three years of operation, we calculate 
the confidence regions for thrice as many supernovae. 
Priors for £l m are Gaussian with an m -prior = 0.05 as be- 
fore, and we also consider <7n m _ pl -i or — 0.015. The lat- 
ter is consistent with the estimated precision of a hypo- 
thetical ground-b ased 1 0° x 10° weak-lensing survey (van 
Waerbeke et al. 1999 ). (As discussed in this reference 



there is a weak dependence of fix in these estimates of 
tt m - We will not pursue this further here.) Uncertainties 
when f2 m is exactly known (elliptic contours) improve 
the expected factor l/v3 as compared with the one-year 
scenario (compare with figure ^). For an £l m prior with 
Cf2 m -prior = 0.05, confidence regions still span consid- 
erable parts of the parameter space. However, with the 
sharper prior, uncertainties in wq and W\ go down to 



Wq 



-1 ± 0.02, Wl 







+0.13 
-0.15 



with an exact A4, and 



are still reasonable when imposing no prior knowledge of 



M: w = -1 ± 0.04, wi = 



+0.15 
-0.17- 



Thus, it seems to us 



that three years of SNAP data backed up with indepen- 
dent high-precision observations of fl m can constrain the 
nature of the dark energy quite well. Note that in the 
above calculations we implicitly assume that the universe 
is known to be flat with high accuracy, since we have im- 
posed fitot = 1. 

Next, we turn our attention to the impact of differ- 
ent redshift distributions of supernovae on the confidence 
region in the (wq,wi) parameter space. The interval z £ 
[0, 2] is divided into 8 subsets containing 500 supernovae 
each uniformly distributed in redshift as: z\ — [0,0.25], 
z 2 = [0.25,0.5], z 3 = [0.5,0.75],.., z 8 = [1.75,2.0]. We 
then compose 4 different experimental situations where in 
each case 2000 supernovae are measured (Am = 0.16 mag 
/SN) sampling events from [zi, Z2, Z3, Za], [zi, Z2, Z7, z$], 
[z 2 , z 4 , z 6 , z 8 ] and [z 5 , z 6 , z 7 , z 8 ], as shown in figure 
the case 9 = (0.3, 0.7, —1, 0) where the mass energy den- 
sity is given a uniform prior with A£l m = 0.1. Clearly, a 
wide range of supernova redshifts is more advantageous 
than only data above or below z = 1. 

4. Effect of adding a small sample of supernovae 

To further illustrate the importance of a small number of 
high-redshift events, we have performed Fisher analyses 
(see appendix A. 2 and (Astier 2001)) to investigate the 
effect of adding 100 supernovae to a large initial sample 
at lower redshift. We do this for initially 2000 supernovae 
uniformly distributed at z £ [0, 1.2]. To emphasize the im- 
portance of events at very low redshift, we do the same ex- 
ercise for initially 2000 supernovae uniformly distributed 
at z £ [0.2,1.2]. Since the effects depend significantly on 
the underlying cosmological model, we investigate thre e 
models: the fiducial model of the SNAP proposal ( 2000 ): 
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Fig. 8. (a) ler and 2a confidenc e regi ons for (u>o,wi) us- 
ing the scenario of Maor et al. ( 2001 ). The elongated el- 
lipses correspond to the assumption of exact knowledge of 
n m , while the larger, non-elliptic regions assume the prior 
knowledge that fi m is confined to the interval f2 m ,truc±0.1. 
Exact knowledge of M. is assumed, (b) 68.3 % and 95 % 
confidence regions for the same cosmology. 



(0.28,0.72,-1,0), a quintessence model derived 



$tru. 

from supergravity considerations (Brax & Martin 1999) 
dtme = (0.2 8, 0.72, -0.8, 0.3), and the model used by Maor 
et al. ( pOOll ) trU c = (0.3,0.7,-0.7,0). 

Figures |lj and |l6| show the effect on the errors of 
fl m and fix when adding 100 supernovae to the samples 
outlined above. As expected, high redshifts pay off when 
determining Vt m and fix, but in case the knowledge of 
M. is poor, it is also important to fill the low-redshift re- 
gion. Note that the curves for exact M. have two minima 
(z max and one intermediate redshift), while those where 
M. is unknown have three (z m i n , Zmax and one interme- 
diate redshift). This is only a manifestation of the fact 
that the optimum redshift distri bution with n parameters 
consists of n 8 functions (Astier 2001 ). (When priors are 
imposed this may no longer be the case.) Furthermore, for 
each curve there are two values of the redshift where it is 
totally ineffectual to add more events. 
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Fig. 9. 68.3 % confidence regions for (wq, W\) in the one- 
year SNAP scenario. The elongated ellipses correspond 
to the assumption of exact knowledge of f2 m : the dash- 
dot-dot-dotted line is with exact M. and the long-dashed 
line corresponds to no knowledge of Ad. The larger, non- 
elliptic regions assume prior knowledge of f2 m : the dash- 
dotted line assumes that fi m is known with a Gaussian 
prior for which t7^ m — prior 

= 0.05; the short-dashed line 
assumes the same prior and exact knowledge of M ; finally, 
the solid line is with f2 m confined to the interval Q m ±0.1 
and exact knowledge of Ai . 
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Fig. 10. 68.3 % confidence regions for (wq, w%) in the one- 
year SNAP scenario without the 100 events for which 
z e [1.2,1.7]. The elongated ellipses correspond to the 
assumption of exact knowledge of Cl m : the dash-dot-dot- 
dotted line is with exact Ai and the long-dashed line cor- 
responds to no knowledge of Ai. The larger, non-elliptic 
regions assume prior knowledge of fl m : the dash-dotted 
line assumes that fl m is known with a Gaussian prior for 
which <7n m -prior = 0.05; the short-dashed line assumes 
the same prior and exact knowledge of Ai; finally, the 
solid line is with Q m confined to the interval fi m ±0.1 and 
exact knowledge of Ai . 
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Fig. 11. 68.3 % confidence regions for (wq, Wi) in the one- 
year SNAP scenario with a constant rate per co-moving 
volume for z <E [0,1.2], and the 100 z E [1.2,1.7] events 
uniformly distributed. The elongated ellipses correspond 
to the assumption of exact knowledge of f2 m : the dash- 
dot-dot-dotted line is with exact M. and the long-dashed 
line corresponds to no knowledge of M. . The larger, non- 
elliptic regions assume prior knowledge of £l m : the dash- 
dotted line assumes that fi m is known with a Gaussian 
prior for which <7Q m — prior 

= 0.05; the short-dashed line 
assumes the same prior and exact knowledge of M. ; finally, 
the solid line is with Q m confined to the interval tt m ±0.1 
and exact knowledge of M. . 
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Fig. 12. 68.3 % confidence regions for (u>o,wi) assum- 
ing the precision of the on e-year SNAP scenario, but the 
cosmology of Maor et al. (2001). The elongated ellipses 
correspond to the assumption of exact knowledge of f2 m : 
the dash-dot-dot-dotted line is with exact M. and the 
long-dashed line corresponds to no knowledge of M. The 
larger, non-elliptic regions assume prior knowledge of f2 m : 
the dash-dotted line assumes that £l m is known with a 
Gaussian prior for which (Jn^— prior 

= 0.05; the short- 
dashed line assumes the same prior and exact knowledge 
of A-f; finally, the solid line is with f2 m confined to the 
interval Q m ±0.1 and exact knowledge of A4. 
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Fig. 13. 68.3 % confidence regions for (u>o,u>i) in the 
three-year SNAP scenario. The elongated ellipses corre- 
spond to the assumption of exact knowledge of f2 m : the 
dash-dot-dot-dotted line is with exact M. and the long- 
dashed line corresponds to no knowledge of A4. The larger, 
non-elliptic regions assume Gaussian prior knowledge of 
f2 m : the dotted line is with ctq m — prior = 

0.05, while the 

dash-dotted line is with crji-prior = 0.015. The solid and 
short-dashed lines assume exact knowledge of M. with the 
same VL m priors as above. 

Figures [l?] - assume a flat universe, and consider 
(wq,Wi) for the same initial distributions. In figures [Tt] 
and [l^ f2 TO is exactly known, while in figures |l^ and 
pfj| a Gaussian prior with o Q m —prior 

= 0.05 is imposed. 
The pay-off with high-redshift events is not as great as 
when determining (£l m ,Qx)- I n particular, note that the 
cosmological-constant model is the worst case of the sce- 
narios we have considered. 



5. Lensing bias 

So far, the analysis has not taken into account any sys- 
tematic errors in the magnitude measurements. However, 
there are several possible mechanisms that can give rise 
to redshift-dependent systematics: attenuation by "gray 
dust" in the intergalactic medium would cause distant 
sources to look fainter than they really are, and evolu- 
tionary effects M = M(z) of the absolute magnitude of 
supernovae la are currently not well-known. Furthermore, 
the effects of gravitational lensing increase with redshift, 
and the corresponding magnitude distributions become 
markedly non-Gaussian for sources at high redshift. 

We have investigated the effects from gravitational 
lensing by using the method of Holz & Wald (199? ), see 
further (Bergstrom et al. 2000). The inhomogeneities are 



modelled as hal os wi th the density profile as proposed by 
Navarro et al. ( 1997). We consider the cosmology exam- 
ined by Maor et al. ( P00l|) , 8 = (0.3,0.7,-0.7,0), and use 
the redshi ft dis tribution given by table 7.2 in the SNAP 
proposal ( 2000 ). Note that this distribution is different 
from the ones used previously. Figure El] shows the lens- 
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Effect of adding 100 SNe to initially 2000 uniformly in (0.2, 1.2) 




Fig. 14. 68.3 % confidence regions for (wq,wi) with 2000 
supernovae. The importance of wide redshift coverage is 
demonstrated by simulating four different synthetic ex- 
periments, all consisting of 2000 SNe: 1234 includes SNe 
uniformly distributed in z G [0,1], experiment 5678 has 
only SNe uniformly distributed in z G [1,2], experiment 
1278 has supernovae in uniformly distributed in two bins, 
z G [0,0.5] and z G [1.5,2]. Finally, experiment 2468 in- 
cludes four bins: z G [0.25, 0.5], z G [0.75, 1], z G [1.25, 1.5] 
and z G [1.75,2]. Clearly, the two experiments with the 
widest redshift coverage provide the best constraints. 
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Fig. 15. The effect on erQ m and an x when 100 supernovae 
are added at a specific redshift z G [0, 2] . The original sam- 
ple consists of 2000 supernovae uniformly distributed over 
z G [0,1.2]. Solid lines correspond to the SNAP fiducial 
model (Q m ,ilx,w ,wi) = (0.28, 0.72, —1, 0), dashed lines 
correspond to (fl m ,fl x ,w ,wx) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (Cl m , fix, Wq, Wi) = 
(0.3,0.7,-0.7,0). 
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Fig. 16. The effect on an m and a^ x when 100 supernovae 
are added at a specific redshift z G [0,2]. The original sam- 
ple consists of 2000 supernovae uniformly distributed over 
z G [0.2, 1.2]. Solid lines correspond to the SNAP fiducial 
model (n m ,Clx,w ,Wi) = (0.28,0.72,-1,0), dashed lines 
correspond to (fl m ,fl x ,w ,wi) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (fl m ,flx,wo,wi) = 
(0.3,0.7,-0.7,0). 

ing effects in the (wo,Wi) space. In this particular pa- 
rameter space the lensing effects are negligible compared 
with the intrinsic uncertainty in the (wo,uii) measure- 
ments. However, sizable effects have to be considered for 
the (fl m , fix) parameter space, especially if fl m contains a 
significant fraction of point-like objects, such as MACHOs 
(Amanullah et al. |200l|). 



6. Discussion 

This analysis stresses the importance of combining inde- 
pendent estimations of the cosmological parameters in or- 
der to probe the nature of the dark energy as accurately 
as possible. For instance, we conclude that a mission for 
observing supernovae over a large redshift range, such as 
the SuperNova/ Acceleration Probe (SNAP), can give rea- 
sonable constraints on the equation of state of the dark 
energy, provided three years of observational data and 
good prior knowledge of the geometry and matter den- 
sity of the universe. To exemplify, we expect SNAP to 
be able to determine the parameters in a linear equation 
of state w(z) — wq + wi z to within ±0.04 for wq and 
t.o.17 f° r w i (one-parameter one-sigma levels), assuming 
a flat universe, the matter energy density known with 
&Q m -prior ±0.015, but no prior knowledge imposed on the 
intercept M.. These estimates assume that the overall er- 
ror budget is not dominated by systematic uncertainties. 
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Fig. 17. The effect on a Wo and a wi when 100 supernovae 
are added at a specific redshift z E [0,2]. The original 
sample consists of 2000 supernovae uniformly distributed 
over z E [0,1.2]. fl m and fix are assumed to be ex- 
actly known. Solid lines correspond to the SNAP fiducial 
model (fl m , fix , wq, wi) = (0.28, 0.72, —1, 0), dashed lines 
correspond to (fl m ,fl x ,w ,wi) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (fl m , fix, wq, w\) = 
(0.3,0.7,-0.7,0). 



Fig. 18. The effect on a wo and a wi when 100 supernovae 
are added at a specific redshift z E [0,2]. The original 
sample consists of 2000 supernovae uniformly distributed 
over z E [0.2,1.2]. fl m and fix are assumed to be ex- 
actly known. Solid lines correspond to the SNAP fiducial 
model (fl m , fix, wo, = (0.28, 0.72, —1, 0), dashed lines 
correspond to (fi m , fl x , w , Wi) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (fl m , fix, wq, wi) = 
(0.3,0.7,-0.7,0). 



With one year of SNAP data, wq could be within 10 % 
provided that the equation of state is assumed to be con- 
stant, W — Wq. 

It is important to realise that data at low as well as 
high redshift is required for an optimal parameter estima- 
tion. Events at very low redshift help to fix the intercept 
Ai , while a wide range of redshifts is needed to break the 
degeneracy in the luminosity distance between different 
cosmologies. 
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Appendix A: Methodology 

We determine two-dimensional confidence regions for sub- 
sets (01,62) € of the parameters 8 = (fl m ,flx,w ,wi), 



while imposing various conditions on the remaining pa- 
rameters. To this end, we construct log-likelihood func- 
tions x 2 based on hypothetical magnitude measurements 
at various redshifts: 



[m(9,M,Zi) - ■m(e tIUC ,M tI n C ,z l )Y 



i=l 



(A.I) 



where m(0,AA, z) is the apparent magnitude of a super- 
nova at redshift z in the cosmology 6 (see section || above), 
and the sum is over bins at different redshifts. The sub- 
script true denotes actual cosmological parameter values. 
The precision Ui of each bin is given by the individual mea- 
surement precision Am and the number of supernovae Hi 
in the bin by cr; = Am/y^. 

Often, we will impose prior knowledge of fl m and/or 
f2 t ot = flm + fix ■ When the parameter of which we have 
prior knowledge is one of the two we are interested in, 
E (0i,02), a Gaussian prior knowledge of with spread 

prior is easily added: 



2 . (0 — fltruc) 
X0 + —2 

8— prior 



(A.2) 



where Xo denotes the x 2 obtained without imposing the 
prior knowledge of 0. In case ^ (0i, 02), we have to inte- 
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Fig. 19. The effect on a Wa and a Wl when 100 super- 
novae are added at a specific redshift z € [0,2]. The 
original sample consists of 2000 supernovae uniformly 
distributed over z 6 [0,1.2]. Ox is assumed to be ex- 
actly known, while m is known within cro m - pr ior = 
0.05. Solid lines correspond to the SNAP fiducial model 
(Q m ,Qx ,w ,wi) — (0.28,0.72,-1,0), dashed lines cor- 
respond to (f2 m , Ox, w ,wi) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (O m , Ox, Wq, Wi) = 
(0.3,0.7,-0.7,0). 



Fig. 20. Effect on a WQ and a Wl when 100 supernovae 
are added at a specific redshift z € [0,2]. The origi- 
nal sample consists of 2000 supernovae uniformly dis- 
tributed over z <E [0.2,1.2]. Ox is assumed to be ex- 
actly known, while TO is known within an m -prior = 
0.05. Solid lines correspond to the SNAP fiducial model 
(O m , Ox, wo, wi) = (0.28,0.72,-1,0), dashed lines cor- 
respond to (O m ,Ox,w ,u>i) = (0.28,0.72,-0.8,0.3), 
and dotted lines correspond to (O m , Ox, wq, w±) = 
(0.3,0.7,-0.7,0). 



grate out 9 from the likelihood L = exp(— ^x ) with some 
prior 7r(6) to obtain Xe-mt- 



Xe-int = -2 In 



J df9exp^-i X 2 j tt(0) 



(A.3) 



Note that the form of (A.3) implies that a constant ad- 
ditive to x 2 simply adds to the integrated log-likelihood 



X0-int' 



-21n 



■W^i»i--(X 2 + A))^(r?) 



Xe-int + A (A.4) 



and that Xe-int — Xe-int min ^ s unaffected by any such con- 
stant. Consequently, we can equally well define 

1 



Xe-int 



-2 In 



'W^Pi-^x'-Xrnin)) 7f(0) 



We will use Gaussian priors 



n(9) = 



2nai 



: exp 



2ai 



-prior 



(A.5) 



(A.6) 



but also uniform priors tt(9) = 1 with 9 confined to an 
interval 9 £ 9 tIUC ± A9. A special case is the treatment of 



the intercept M, for which we assume both exact knowl- 
edge, but also no prior knowledge at all. Hence, integrating 
M. over all possible values M. E (—00,00), we ob tain an 
analytic expression for XM-int> see appendix |A.l[ 

Given the appropriate x 2 function, 68.3 % and 95 % 
confidence regions are defined by the conventional two- 
parameter x 2 levels 2.30 and 5.99, respectively. Similarly, 
one-parameter one- and two-sigma levels correspond to 
X 2 = 1 and 4, respectively. In some cases we need to cal- 
culate x 2 f° r three parameters, and subsequently project 
onto the (9i,9 2 ) plane of interest. This can be done by 
setting x 2 = min [x 2 (' ■ ■ , $3)], where the minimisation of 
X 2 is performed with respect to variation of #3. Confidence 
regions for (9i, O2) can then be determined using the usual 
two-parameter x 2 levels. 



A.l. Integration over the intercept M. 

When the intercept is assumed to be exactly known A4 = 
-Mtrue, it will cancel in the expression for x 2 , so that we 
obtain the log-likelihood x 2 as 



E 



A 2 



(A.7) 
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A. 2. Fisher matrix analysis 
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Fig. 21. 68.3 % confidence regio n for (wp,wi) for one 
year of SNAP for the cosmology in |2001 , central values at 
Wo = —0.7 and w± = (solid line). The dashed curve in- 
corporates the bias effects form gravitational lensing mag- 
nification. 



A = 51og 10 [d' L (0,Zi)] - 51og 10 [d' L (e tT ue, Zi)] 



(A. 



Note that x 



0-2 



x 2 (0t 



by construction. 



If no prior knowledge of Ai at all is assumed, we can in- 
tegrate the general \ 2 function (A.l) over M. 6 (— oo,oo) 



to obtain an analytic expression for \ = xjU-i 

1 



X z = -2 In 



B 



dMexp ( -^x 2 



C 



C 
- A 

L—l 1 



2tt 



C 



n 

Ei- 



(A.9) 
(A.10) 
(All) 

(A.12) 



Note that this expression is independent of -Mtrue> and 
that we imposed a uniform prior tt(M.) = 1 in the inte- 
gration. It is also worth pointing out that 



Xmin = X (Otvue) = In 

More importantly, 



x 2 



Xmin 



x 



c 



<x 2 



X 



Xmin' 



(A.13) 



(A.14) 



where the equality holds when B = 0. Note that this is 
the case not only when 9 — 8 truc , but in general also on 



a hypersurface in parameter space. The inequality (A.14) 
ensures the intuitive notion that X 2- Xmin contours always 



should lie outside corresponding \ 



c2 



contours. 



For efficient estimators (i.e., in the large sample limit), we 
can obtain the Fisher matrix by finite-difference evalua- 
tion of the expression 



F. 



d 2 log(L) 



dOjdOk 



i d 2 x 2 



2 dOjdOk 



(A.15) 



where, with negligible bias, we can take 6 = ©true- The 
covariance matrix is now given by the inverse of F. 

In the quadratic approximation of \ 2 (with \ 2 based 
on the luminosity distance rather than the apparent 
magnitude m), the Fisher matrix is obtained as 



F jk = J2 h j( Z i) h I( Z i)> 



hj(zi) 



1 ^± 

ui 89 j 



(A.16) 



(A.17) 



where the precision can be expressed in terms of the rela- 
tive precision p as cr, = pd^Zi). It is straight-forward to 
add prior knowledge of any combination of the parame- 
ters 6. Imposing no prior knowledge of M. corresponds to 
letting the scale of d^ be unknown: d^ = Q d' L . 



It should be noted that, even though equation (A.16) 
is an approximation, it gives uncertainties in accordance 
with the analysis in section p| (co mpare, for instance, max- 
imum values in figures [l^ - |20| with relevant cases in ta- 
bles Al and |A.3j). In addition, for inefficient estimators 



(i.e., non-ellipsoidal confidence regions), the approximate 
Fisher analysis roughly gives the mean errors of parame- 
ters. 
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±0.021 


±0.053 


( ^ m j 


fix) (constant rate/volume at 2 £ [0, 1.2]) 


exact Al, no prior fl m 

exact Al, Gaussian fl m , <Jn m - pT i OI — 


0.05 


±0.016 
±0.015 


±0.030 
±0.028 






no prior Al, no prior fl m 

no prior Al, Gaussian fl m , an^ -prior 


= 0.05 


±0.017 
±0.016 


+0.079 
-0.087 
+0.077 
-0.084 



Table A.l. One-parameter one-sigma ranges for (Sl m ,Slx) in the one-year SNAP scenario. The quoted parameter 
ranges for a parameter 9 are obtained by finding the extremal values of 9 for which \ 2 = 1 . 



(fi m , wo), fixed fl tot = 1 exact Al, no prior fl m -0M9 ±0.048 

exact Al, Gaussian fl m , an m - P rior = 0.05 to'.oi? ±0.045 

no prior Al, no prior fl m to.023 ±0.078 

no prior Al, Gaussian fl m , an m -prior = 0.05 i'um _ ±0.071 

(fl m ,w ), fixed fix = 0.72 exact Al, no prior fi m -"i.017 ~ ±0.030 

exact Al, Gaussian fi m , <7n m - prior = 0.05 — o'.oii ±0.028 

no prior Al, no prior fl m to.021 ±0.055 

no prior M, Gaussian fi m , o"n m -prior = 0.05 ilhoi" ±0.052 

(fix, w ), fixed fl m = 0.28 exact Al, no prior fix - +q 010 ±0.12 

exact Al, Gaussian fix, otj x -prior = 0.05 -0.010 ±0.045 

no prior Al, no prior fix -0.026 ±0.12 

no prior Af, Gaussian fix, crn x -prior = 0.05 - -0.024 ±0.049 

(fix, wo), Gaussian fl m, C"Q m — prior = 0.05 and Gaussian fi to t, & Otot —prior = 0.05, exact Af ±0.07 

(fix, wo), Gaussian fi m , crn m - pr ior = 0.05 and Gaussian fitot , o"n tot - pr ior = 0.05, no prior Af - -o'.ot ±0.09 



Table A. 2. One-parameter one-sigma ranges for (Q m ,wo) or (Qx,wo) in the one-year SNAP scenario. The quoted 
parameter ranges for a parameter 9 are obtained by finding the extremal values of 9 for which \ 2 — 1, with the 
additional requirement w > —1. 
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(WQ,Wl), fitot = 1 


exact M, exact f2 m 


±0.031 


±0.14 


exact M, fi,n £ Om.truc ± 0.1 


+0.13 
-0.066 


+0.48 
-0.76 




exact M, Gaussian fl m , dfi m _ pr i or = 0.05 


+0.065 
-0.052 


+0.31 
-0.46 




no prior M, exact fl m 


±0.064 


±0.18 




no prior M, Gaussian fi m , <7f2 m - P rior = 0.05 


+0.077 
-0.074 


+0.35 
-0.53 



(wo,wi), (no z £ [1.2, 1.7] events) 


exact Ai, exact fl m 


±0.034 


±0.16 




exact M, Q.m £ fim,truc ± 0.1 


+0.13 
-0.085 


+0.50 
-1.03 




exact M, Gaussian S7 m , <Jn m - pTioT — 0.05 


+0.068 
-0.059 


+0.32 
-0.48 




no prior M, exact O m 


±0.070 


±0.21 




no prior M, Gaussian Q m , aci m - pT i or — 0.05 


+0.085 
-0.083 


+0.36 
-0.54 



(wq,Wi), (constant rate/volume at z £ [0, 1.2]) 



exact M, exact Q, m 

exact M, Q.m £ fim.truc ±0.1 

exact M, Gaussian Q, m, ffl m - prior = 0.05 

no prior M, exact f2 m 

no prior M, Gaussian Q. m , o-f2 m -prior = 0.05 



±0.038 

+0.13 
-0.066 
+0.065 
-0.054 
+0.14 
-0.15 
+0.14 
-0.15 



±0.16 

+0.50 
-0.96 
+0.33 
-0.53 
+0.27 
-0.26 
+0.42 
-0.66 



(wo,wi) (Maor et al. cosmology) 



exact M, exact f2 ra 

exact M, Q.m £ fim.truc ±0.1 

exact M, Gaussian Q m , dfi m _ pr i or = 0.05 
no prior M, exact fi m 

no prior M, Gaussian Q m , <7f2 m -prior = 0.05 



±0.028 

+0.10 
-0.10 
+0.054 
-0.052 

±0.057 

+0.070 
-0.072 



(wo,wi), fitot = 1j (three-year SNAP) 



exact M, exact fi m 
exact M, Gaussian f2 
exact M, Gaussian Q, m , usi„ 
no prior M, exact f2 m 
no prior M, Gaussian Q m , ffsi„ 
no prior M, Gaussian Q m , aci^ 



-prior 
-prior 



0.05 
0.015 

■ = 0.05 
. = 0.015 



±0.018 

+0.060 
-0.038 
+0.023 
-0.024 

±0.036 

+0.062 
-0.047 
+0.038 
-0.039 



±0.11 

+0.26 
-0.63 
+0.16 
-0.22 

±0.16 

+0.20 
-0.25 



±0.081 

+0.29 
-0.36 
+0.13 
-0.15 
±0.11 
+0.31 
-0.42 
+0.15 
-0.17 



(wo,wi) (Maor et al. scenario) 



exact M, exact fi m ±0.014 ±0.044 

exact MO (= O +0 1 +0.094 +0.21 

CAd ^ 1 ' l ; iL m c itm,truc -J- W.J- —0.051 —0.34 

Table A. 3. One-parameter one-sigma ranges for (wq,wi) in the one-year SNAP scenario. Note that the two last 
sections instead refer to the three- year SNAP scenario and the scenario of Maor et al. ( [2001 ), respectively, also 



discussed in section B.3. The quoted parameter ranges for a parameter 9 are obtained by finding the extremal values 
of 6 for which x 2 = 1. 



